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We theoretically study the spin fluctuation and superconductivity in Lallll and Smllll iron-based 
superconductors for a wide range of electron doping. When we take into account the band structure variation 
by electron doping, the hole Fermi surface originating from the d X 2_ y 2 orbital turns out to be robust against 
electron doping, and this gives rise to large spin fluctuations and consequently s± pairing even in the heavily 
doped regime. The stable hole Fermi surface is larger for Smllll than for Lallll, which can be considered 
as the origin of the apparent difference in the phase diagram. 

KEYWORDS: iron pnictide superconductors, Fermi surface, spin fluctuation, first-principles band calculation 



The pairing mechanism of the iron-based superconduc- 
tors has been of great interest ever since its discovery. 4 ) 
The possibility of spin-fluctuation-mediated pairing in 
the iron-based superconductors has been proposed from 
the very beginning of the study. 2,3 - 1 The key point here is 
the presence of the disconnected electron and hole Fermi 
surfaces, which gives rise to spin fluctuations around the 
wave vector that connects the Fermi surfaces. The result- 
ing pairing state is the so-called s± state in which the 
sign of the gap function changes between electron and 
hole Fermi surfaces. 

Recent observation of superconductivity in "1111" 
compounds LnFeAsOi-xR x (in=La,Sm, etc.) for a wide 
range of x ' 5 ' has raised an interesting issue regarding the 
possibility of spin-fluctuation-mediated pairing. There, 
the superconducting state persists up to a large x of 
~ 0.4, for which hole Fermi surfaces are expected to be 
completely wiped out in a rigid band picture. Therefore, 
whether the spin fluctuation can be responsible for pair- 
ing even in such a heavily doped regime is certainly an 
intriguing problem. 

To address this issue, in the present paper we study 
the spin fluctuations and superconductivity of Lallll 
and Smllll by applying random phase approximation 
(RPA) to a five orbital model that takes into account 
the effect of the variation of the band structure due to 
electron doping. It is found that the hole Fermi surface 
around the wave vector (jr, n) that originates from the 
dx^-Y 2 orbital is robust against electron doping, and 
this induces large spin fluctuations and s± pairing even 
in the heavily electron doped regime. The stable hole 
Fermi surface is larger for Smllll than for Lallll, and 
we propose that this is the origin of the difference in 
the experimentally observed T c vs. x phase diagram, in 



which the latter material has a two-dome feature, while 
the former has a single-dome shape. 4, 5 ) 

We start with the first principles band calculation for 
the two materials. The band calculation is performed us- 
ing the VASP package 6 -* for Lallll and Smllll com- 
pounds. Since the rigid band picture is expected to be 
unreliable especially in the heavily doped regime, we take 
into account the effect of the band structure variation due 
to doping by (i) adopting the lattice structure parame- 
ters determined experimentally for each doping rate x, 
and (ii) adopting virtual crystal approximation in which 
the oxygen and the fluorine pseudopotentials are mixed, 
assuming that the hydrogen doping has an affect of in- 
creasing the average valence of the oxygen site by +.T. 4 -* 
We construct a five orbital tight-binding Hamiltonian in 
the unfolded Brillouin zone 3 ) exploiting the maximally 
localized Wannier functions. 7 - ) The five Wannier orbitals 
have five different symmetries (dxy-, dyz, dzx, d 3 z 2 -n 2 
and dx 2 -Y 2 )i where X, Y refer to the direction rotated 
by 45 degrees from the Fe-Fe direction x, y. 

We construct a many body Hamiltonian by taking into 
account the multi-orbital electron-electron interaction, 
and apply random phase approximation (RPA). 9 ' 10 ) Us- 
ing the Green's function Gi m (k) (fc = (fe, iu n )), which 
is a 5 x 5 matrix in the orbital representation, the irre- 
ducible susceptibility matrix is given as 



xl MMi {q)=Y J G llh {k + q)G ul2 {k) 1 



k 



(1) 



and the spin and the charge (orbital) susceptibility ma- 
trices are obtained as 



x*(q) = 



x 



l-Sx°(q)' 



(2) 
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Fig. 1. The Fermi surface for various doping rates for Lallll 
(top) and Smllll (bottom). The arrows 1 and 2 indicate the in- 
teractions between the Fermi surfaces which give rise to spin fluc- 
tuations and superconductivity. 
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Fig. 2. The Stoner factor for Lallll plotted against the doping 
rate. 
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where S and C are the corresponding interacting vertex 
matrices. These matrices have l\ to I4 (k — 1,...,5) as 
orbital indices. We calculate the Stoner factor defined as 
the largest eigenvalue of the matrix §x°(q)- This qual- 
itatively measures the strength of the low energy spin 
fluctuations, and a tendency toward magnetism. 

The Green's function and the effective singlet pairing 
interaction, 

V S (q) = ls X s(q)S - \cUq)C + ^(S + C), (4) 

are plugged into the linearized Eliashberg equation for 
superconductivity, 



MhiA k ) 



N a^ z^ v hhhh(q) 



x Ghh ( k ~ q)(t>hh (k - q)G hk (q - fc)(5) 
The 5x5 matrix gap function cf>i m in the orbital repre- 



sentation along with the associated eigenvalue A is ob- 
tained by solving this equation. Since A(T) = 1 signals 
T = T c , the eigenvalue calculated at a fixed temperature 
is a qualitative measure for T c . We present the gap func- 
tion transformed into the band representation with a uni- 
tary transformation. The RPA calculation is performed 
for 64 x 64 x 4 sites and 2048 Matsubara frequencies, and 
the temperature is T = 0.02eV. The electron-electron in- 
teraction is taken from ref.8, but we multiply all of them 
by a factor of / = 0.42 since RPA overestimates the ten- 
dency towards magnetism. 

We first show in Fig.l the Fermi surface for various 
values of the doping rate x. The hole Fermi surfaces (ai, 
012) around the wave vector (0, 0) originating from the 
dxz/dyz orbitals shrink, while the electron Fermi sur- 
faces (/3i, (3 2 ) around (0, 7r)/(7r,0) become larger with 
doping, as expected. What is interesting is that the vol- 
ume of the hole Fermi surface (7) around (re, re), which 
originates from the d X 2^y 2 orbital, barely changes with 
doping for both materials. This is due to the band struc- 
ture variation with doping, whose origin can be two 
folded. One is the lattice structure variation, and the 
other is the change in the electric charge of the FeAs and 
LnO layers. Both effects are taken into account in the 
present band calculation, and it turns out that the lat- 
ter effect actually dominates. Namely, the energy level 
of dxz/dyz orbitals (that are oriented toward the LnO 
layers) is pushed down relative to c?x 2 -y 2 owing to the 
increase of the positive charge in the LnO layers. Com- 
paring the two materials, the 7 Fermi surface is larger 
for Smllll due to higher pnictogen position. 11 ' 

In Fig. 2, we show the Stoner factor as a function of x. 
It can be seen that the Stoner factor actually increases 
with doping for a fixed electron-electron interaction, as 
opposed to a naive expectation based on "rigid band + 
Fermi surface nesting" picture. To see this in more detail, 
we show in Fig. 3 the ratio of the "intra-orbital spin sus- 
ceptibility" between dx^-yz and dxz/dyz orbitals. 11 ' 
Here, the intra-orbital spin susceptibility is defined as the 
diagonal element of the spin susceptibility matrix that 
has the same orbital indices. In Lallll, the dxz/dyz 
susceptibility dominates for small doping, but as the dop- 
ing rate is increased, dx^-Y? tends to dominate, par- 
ticularly around the wave vector (re, 7r/3)/(7r/3,7r). 13 ' 
The contribution to this wave vector comes from two 
interactions. One from the 7-/3 (electron-hole, arrow 1 
in Fig.l) interaction, and another from /3i-/?2 (electron- 
electron, arrow 2) interaction. 11 ' 14 ~ 16 ' These two contri- 
butions, both coming mainly from the intra-orbital re- 
pulsion within the dx 2 -y 2 orbitals, cooperate to give a 
large enhancement in the incommensurate spin fluctua- 
tions. In fact, the enhancement of the spin fluctuation 
around this wave vector is consistent with a recent ex- 
perimental study. 12 ' The cooperation occurs because the 
two contributions accidentally have similar incommen- 
surate wave vectors especially for large a;(~ 0.3 — 0.4). 
As for the material dependence, the difference between 
Lallll and Smllll can be seen by comparing the up- 
per and lower panels in Fig3. While a crossover from 
dxz/dyz to dx 2 -Y 2 dominating regimes is seen in La, 
the dx 2 -Y 2 contribution is large through the entire dop- 
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Fig. 3. The ratio of the intra-orbital spin susceptibilities X s X 2_ V 2/Xxz/YZ ^ or var i° us doping rate for Lallll(top) and 
Smllll(bottom). 



ing regime for Smllll. This is due to the presence of the 
larger 7 Fermi surface in the latter material. 

The two origins for the spin fluctuations, i.e., the 7-/? 
and /3i-/32 interactions cooperate toward antiferromag- 
netic instability, but they have competing effects regard- 
ing the occurrence of superconductivity. Namely, as dis- 
cussed in previous studies, 3, 11,14 ~ 16 ) the 7-/3 interaction 
leads to s± pairing, while the (3\-(32 interaction favors 
(i-wave pairing because this interaction acts to change 
the gap sign between /3i and @2- To see which one of 
the interactions dominate in the pairing interaction, we 
calculate the eigenvalue of the Eliashberg equation for 
s and d-wave pairings. As shown in Fig. 4 (a), s-wave al- 
ways dominates over <i-wave, and the gap function has 
the fully gapped s± form as shown in Fig. 4(b), where the 
gap changes sign between the electron and hole Fermi 
surfaces. 

The present result suggests that the dx 2 -Y 2 hole 
Fermi surface plays an important role in the occurrence 
of superconductivity even in the heavily doped regime. 
One should note, however, that the nesting (in the con- 
ventional sense of the term) between the electron and the 
hole Fermi surfaces is not good, since the volume is very 
different. Nevertheless, the pairing is strongly enhanced. 
This is because the portions of the bands away from the 
Fermi level contribute to the pairing interaction through 
finite energy virtual processes. Therefore, one should be 
careful about using the word "Fermi surface nesting" as 
the origin of the pairing interaction, at least in the naive 
sense of the term. 

The material dependence seen in the intra-orbital con- 
tribution to the spin susceptibility is also reflected to su- 
perconductivity. In Fig. 4(c), we plot the eigenvalue of the 
Eliashberg equation against the doping rate for the two 
materials, normalized by its value at the lowest doping 
rate. We see that it is strongly enhanced in the heav- 
ily doped regime for Lallll, while the variation is rela- 



tively weak for Smllll. This is because the contribution 
weight to the spin fluctuations from the dx 2 -Y 2 orbital 
is less dependent on the doping rate in Sm compared to 
that in La. The weak doping dependence of A in Smllll 
is qualitatively consistent with the small variance of T c 
for a wide range of x observed in Smllll. 5 ) Although 
a complete understanding of the two-dome shape of the 
Lallll phase diagram lies beyond the present theoreti- 
cal approach (which does not take into account the self 
energy and the vertex corrections), smaller eigenvalues 
obtained for small x in Lallll suggests a fragile nature 
of the superconductivity in this regime. 

To summarize, we have analyzed the spin fluctuation 
and superconductivity of the 1111 iron pnictides in a 
wide range of electron doping rate by taking into account 
the band structure variation within the virtual crystal 
approximation. The dx 2 -Y 2 hole Fermi surface around 
(71", 7r) is found to be robust for the entire doping regime, 
and plays an important role on both the spin fluctua- 
tions and superconductivity. In Lallll, the dxz/dyz 
orbitals have relatively large contribution in the small x 
regime, while in Smllll, the dx 2 -Y 2 dominates over the 
entire regime. This difference in the dx 2 -Y 2 contribution 
is plausibly the origin of the difference in the phase di- 
agram between the two materials. After completion of 
this study, a paper that studies a similar problem for 
LaFeAsOi-^Hz has been posted on the arXiv. 17 ) 
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Fig. 4. (a) The eigenvalue of the Eliashberg equation as a func- 
tion of the doping for Lallll. (b) The gap function for the s-wave 
pairing for band 3 (with the hole Fermi surfaces) and band 4 (the 
electron Fermi surfaces). The solid lines represent the Fermi sur- 
faces, (c) The normalized eigenvalue for the s-wave pairing against 
the doping rate for La and Sm. The values are normalized by the 
values at x = 0.13 to focus on the x dependence. 
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